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Abstract. A theoretical analysis of the lateral displacement (Goos-Hanchen shift) 
of spatially limited beams reflected from Attenuated Total Reflection (ATR) devices 
in the Otto configuration is presented when backward surface plasmon polaritons are 
excited at the interface between a positive refractive index slab and a semi-infinite 
metamaterial with negative refractive index. First, the stationary phase approximation 
and a phenomenological model based on the properties of the complex poles and zeroes 
of the reflection coefficient are used to demonstrate that: i) the excitation of backward 
surface waves can lead to both negative and positive (and not exclusively negative) 
Goos-Hanchen shifts, and ii) the sign of the shift depends on whether the value of 
the coupling layer thickness is higher or lower than a critical value characteristic of 
the ATR structure. Second, these findings are verified through rigorous calculations 
of the spatial structure of the reflected beam. For incident beams with a Gaussian 
profile, the lateral shift calculated as the first moment of the field distribution of the 
reflected beam agrees quite well with the predictions of approximate analysis. Near 
the resonant excitation of the backward surface plasmon polariton, large (negative 
or positive) Goos-Hanchen shifts are obtained, along with a splitting of the reflected 
beam. 


RAGS numbers: 81.05.Xj,73.20.Mf,78.68.-bm 

1. Introduction 

The Goos-Hanchen (GH) lateral shift is an optical phenomenon in which a spatially 
limited beam reflected at an interface near conditions of total internal reflection is 
displaced from the incident beam position. Althongh conjectnred by Isaac Newton 
in the 18th centnry [1], experimental evidence was only provided in 1943 by Hilda 
Hanchen (a doctoral stndent at the time) and her thesis advisor Fritz Goos |2l |3] . The 
effect was theoretically explained by Artmann [1], who noted that each plane wave 
component in the expansion of the incident beam undergoes a slightly different phase 
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change after total internal reflection, so the sum of all the reflected components produces 
a lateral displacement of the reflected beam from the geometrical optics prediction. 
When the stationary phase approximation is used to calculate the sum of all the reflected 
components, the lateral beam shift is given by 
dip 


A = - 


da' 


( 1 ) 


where a is the tangential component of the wave vector in the medium of incidence and 
ip is the phase difference between the reflected and incident waves. Incidentally, we note 
that Artmann’s argument shows that a prerequisite for the existence of a Goos-Hanchen 
shift is not the excitation of surface or evanescent waves but rather the phases of the 
different plane wave components in the expansion of the incident beam being shifted 
differently upon reflection. 

The GH lateral shift and other nonspecular phenomena [5], such as beam prohle 
deformation, focal displacement, spatial waist modihcation and angular deviation, have 
been theoretically investigated in diverse conhgurations which include multilayers [6] and 
periodic structures [7]. For beams reflected from a single flat interface, the GH effect 
is much smaller than for beams reflected from periodically corrugated or multilayered 
structures, in which the shift could be of the same order of magnitude as the beam width 
[8]. Apart from the case of linear isotropic dielectric media, nonspecular phenomena 
have been estimated and measured for other kind of homogeneous materials, such as 
chiral |9], anisotropic uniaxial [inuniiis] or dielectric-magnetic media with negative 
refraction index da [HE]. A unihed linear algebra approach to study the generalized 
shifts for the polarization components of reflected light beams is presented in dB]- For 
a self-consistent description of the GH shift and related beam-shift phenomena together 
with an overview of their most important extensions and generalizations, the reader is 
referred to im. 

For total internal reflection from a boundary between two positive index media, the 
GH shift is typically the same order of magnitude as the wavelength and positive (i.e., 
the reflected beam is shifted to the other side of the normal from the incident beam). 
On the contrary, for reflection from a metallic boundary and provided that the beam 
is p polarized, the GH shift may be negative (i.e., the reflected beam can be shifted 
towards the same side of the normal as the incident beam). As shown in Ref. [18], this 
negative shift is due to the fact that the p polarized evanescent wave in the metal is 
backward, i.e., it has an energy flow opposite to the direction of the phase velocity along 
the interface. For total internal reflection at the boundary between a positive index 
medium and a negative index medium [131 El E ]5 lateral shift is also negative 
because the evanescent wave is again backward [T9] . 

Although the beam shift at a single flat interface can be either positive or negative 
depending on the type of the surface waves (forward or backward) excited by the 
incoming beam, the same is not true for beam reflection from other structures. For 
example, it is well known (see [20] and references therein) that the excitation of a forward 
surface wave in a prism-waveguide coupling system can result in either a positive or a 
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negative GH shift. However, the erroneous belief that the excitation of a backward 
surface wave always results in a negative GH shift sometimes appears in the literature, 
particularly when multilayers with negative index materials are involved [211 ESI EH]. 

The purpose of this paper is twofold: (a) to study the lateral displacement of 
spatially limited beams reflected from Attenuated Total Reflection (ATR) structures 
under conditions of resonant excitation of backward surface plasmon polaritons (SPPs) 
and (b) to apply the phenomenological model presented in [ 2 l] to clearly evidence 
that the excitation of backward surface waves in multilayers with negative index 
metamaterials can lead to both negative and positive (and not exclusively negative) 
GH shifts, depending on whether the thickness of the coupling layer is higher or lower 
than a critical value representing a real zero of the reflection coefficient. 

The plan of this paper is as follows. In Section [2] we present some technical 
background and two alternative methods that will be used to numerically calculate the 
GH lateral shift, namely the stationary phase approximation and the rigorous calculation 
of the held distribution of the rehected beam. In Section |3] we discuss the GH shift near 
conditions of resonant coupling between the incident held and backward surface waves 
in terms of the behavior of the complex poles and zeroes of the rehection coefficient 
[23]. We give examples that show the existence of a critical thickness of the coupling 
layer: above this critical value the lateral shift is positive, whereas below this critical 
value the lateral shift is negative. In Section ID we present rigorous calculations of 
the transverse distribution of the rehected electric held near the resonant excitation of 
the backward surface plasmon polariton. We investigate the inhuence of the angle of 
incidence and the thickness of the ATR device for the optimal condition of coupling 
between the incident radiation and backward SPPs. Large (negative or positive) GH 
shifts are obtained, accompanied by a splitting of the rehected beam. In Section [5] we 
summarize and discuss the obtained results. 

2. GH shift calculation 

Since the spatial periodicity associated with surface plasmon polaritons is less than the 
spatial periodicity which could be induced by an incident photon on the surface, surface 
plasmon polaritons cannot be resonantly excited by plane waves. This difficulty can be 
overcome by using ATR techniques, which require the introduction of a second surface, 
usually the base of an isosceles prism. In the Kretschmann conhguration, the prism 
is positioned against the plasmonic medium (metal o metamaterial), while in the Otto 
conhguration [25] the prism is positioned very close to the surface of the plasmonic 
medium. Figure [H shows an Otto-ATR structure illuminated by a spatially limited 
gaussian beam (width 2a). The beam is incident upon the base of the prism a.t y = 0, 
the center of the beam waist is located at a: = 0 and the angle of incidence of the 
beam axis is Oq. Medium 3 is a magnetic metamaterial with complex values of electric 
permittivity £3 and magnetic permeability These constitutive parameters satisfy the 
condition for power how and phase velocity in opposite directions [26]. Medium 1 (the 
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Figure 1. Schematic diagram of the system. It shows the lateral shift A of a bounded 
beam at an interface. The dashed lines show the boundaries of the reflected beam as 
predicted by geometrical optics theory. 


prism) and medium 2 (the coupling layer, thickness d) are nonmagnetic (pi = /i 2 = 1 ) 
and have real and positive electric permittivities (^i and £2 respectively). When the 
values of the relative constitutive parameters es /£2 and /i 3 //i 2 are properly chosen |27] . 
the interface 2-3 between the coupling layer and the metamaterial can support SPPs 
with time-averaged Poynting vector directed opposite to the phase velocity (backward 
SPPs). The examples of this paper are for the case (£ 3 / 62 ) < —1, that corresponds 
to s-polarized backward SPPs. 

When s-polarized incident radiation reaches the base of the prism with an angle 
6q greater than the critical angle of total reflection, the evanescent held can resonantly 
couple with the s-polarized backward SPP of the interface 2-3. We assume that the 
electric held of the incident beam is given by 

Ei = Ei{x,y = 0) z = exp [-(—)^ ■PiaQx]z, (2) 

with ao = sin 6 * 0 , c the speed of light in vacuum, and Ox = a sec( 6 'o). An 

exp {—iut) time dependence is implicit, with u the angular frequency, t the time and 
i = \/—1. Using a Fourier integral representation, the incident held can be written as 

/ OO 

A{a) exp {i ax) da, (3) 

-OO 

where 


A{a) = 0=exp[-{a-ao)H^)^], (4) 

is the angular spectral distribution of the incident beam. The electric held of the 
rehected beam at the interface y = 0 is given by 

/ OO 

R{a) A{a) exp {i a x) da (5) 

-OO 

where R{a) = |i?(a)| expiip{a) is the complex rehection coefficient of the multilayer as 
a function of the spectral variable a = sin 6. The beam shift can be obtained 
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by finding the location where \Er{x,y = 0)| is maximal. Using the stationary-phase 
approximation and assnming a linear variation for (p{a) and that the beam experiences 
total internal reflection, this location is given by eqnation ([T]). However, if ip is not a 
linear fnnction of a or if |i?(a)| 7 ^ 1 along the spectral width of the incident beam, the 
reflected beam may snffer severe distortions and the stationary-phase result given by 
equation ([1]) is usually unsuitable to describe the lateral shift correctly [28]. In such 
cases, it is preferable to define the GH shift using the normalized first moment of the 
electric field of the reflected beam m 


. ( 1 ) _ IZc^\^r(x,y = 0)\‘^dx 
iZc \Er{x,y = 0)\‘^dx 


( 6 ) 


In this definition, the shift is expressed as an average of the x-position of the reflected 
beam on the surface of interest and no approximation is done over the amplitude and 
phase of the reflected plane wave components of the beam. When the spectral width 
of the incident beam is small, it is possible to derive equation © as a first order 
approximation of equation O I3D1. 


3. Zero-pole model for the GH shift 

To numerically illustrate the lateral displacement of spatially limited beams under 
conditions of resonant excitation of backward SPPs, we choose an Otto-ATR structure 
with El = 2.25, /ii = 1, £2 = 1, P 2 = 1, £3 = —1.6 -|- 0.001 i and p .3 = —0.7 -|- 0.001 i 
For these constitutive parameters, the interface 2-3 between the coupling layer and the 
metamaterial supports s-polarized SPPs with time-averaged Poynting vector directed 
opposite to the phase velocity [27] . 



Re z 

Figure 2. (Color online) Trajectories of the zero (zo) and the pole {zp) in the complex 
plane as functions of d/X for an Otto - ATR structure with ei = 2.25, yi = 1, £2 = 1, 
^2 = 1, £3 = -1.6-1-0.001 i and ys = -0.7-1-0.001 i. For these constitutive parameters, 
the interface 2-3 supports s-polarized, backward SPPs. 
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As shown by equation ([S]) the distortions of the reflected beam are ruled by the 
behavior of R{a) = |-R(q!)| expi<y9(Q!). The exact form of this reflection coefficient is 
given for instance by equation (1) in [2l]. However, to clearly evidence the facts that 
i) the excitation of backward surface waves can lead to both negative and positive (and 
not exclusively negative) GH shifts and that ii) the sign of the GH shift depends on 
wheater the value of the coupling layer thickness is higher or lower than the critical 
value characteristic of the ATR structure, hrst we choose to represent this quantity 
within the frame of a phenomenological model based on the properties of the complex 
poles and zeroes of the reflection coefficient. According to this model, used previously 
for the case of Kretschmann-ATR structures with metamaterials m, the reflection 
coefficient R can be written as 

= (T, 

where z is the analytic continuation of sin^ to the complex plane, zo and Zp denote 
respectively the complex zero and the complex pole of R, A is the wavelength in 
vacuum and ({z, d/\) is a complex regular function near zq and Zp that does not change 
signihcantly near Zp. 



e (deg) 


Figure 3. Phase y of the reflection coefficient for an Otto-ATR structure with 
Si = 2.25, = 1; £2 = 1; = 1; £3 = —1.6+ 0.001 1 , ^3 = —0.7+ 0.001 f and 

d/X = 0.7. 

As explained in [23] (see hgures 2 and 3 therein), one of the main characteristics of 
using the zero-pole model to represent the optical response near the excitation of SPPs 
is that it provides a direct visualization of the behavior of the phase curves (p{a) in terms 
of the spectral variable a. This characteristic makes the zero-pole model particularly 
suited for the study of the lateral GH shift, since the quantity that hxes the sign of the 
GH shift in the stationary phase approximation equation ([T|) is precisely the derivative 
of the phase. 
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The interested reader is referred to [23] for further details, here we limit ourselves 
to mentioning the key result regarding the sign of the GH shift, namely that phase 
curves (p{a) have a very different behavior depending on the location in the complex 
plane of both the zero and the pole: when < 0, the phase curve (p(a) is a 

monotonically increasing function, whereas when 7s{zp) > 0, the phase curve is a 

monotonically decreasing function, exhibiting hrst a maximum and then a minimum. 
In terms of the sign of the GH shift, negative GH shifts should be expected in the hrst 
case, whereas positive GH shifts should be expected in the second case. 



0 (deg) 


Figure 4. Phase ip of the reflection coefficient for an Otto-ATR structure with 
El = 2.25, Pi = 1, £2 = 1, fJ ,2 = 1) £3 = —1-6+ 0.001*, ^3 = —0.7 + 0.001* and 
d/X = 0.9. 

Figure 121 where we plot the parametric trajectories of zo{d/\) and Zp{d/X) 
calculated for the values of the constitutive parameters used in this example, clearly 
shows that both cases (Q'(^o) < 0, negative GH shift, or Q'(zo) ^(^^p) > 0, positive 

GH shift) can occur near the excitation of a backward SPP. The pole trajectory cannot 
cross the real axis, (if it did, inhnite rehectance would result for a real angle of incidence), 
but the trajectory of the zero of the rehection coefficient is not limited and in this 
example it crosses the real axis for a critical value oi d/X ~ 0.8. The value of Zq 
0.74) at this critical thickness of the coupling layer corresponds to an angle of 
incidence 60 = 47.8° for which a plane wave is totally absorbed by the ATR structure. 
We conclude that, for values of d/X lower than this critical value, the zero-pole model 
predicts that the function (p{a) is monotonically increasing (as the curve shown in hgure 
[3] calculated for d/X = 0.7), and according to the stationary phase approximation, 
negative GH shifts are to be expected. On the other hand, for values of d/X above 
this critical value, the zero-pole model predicts that the function (p{a) is monotonically 
decreasing (as the curve shown in hgure 01 calculated for d/X = 0.9), and according to 
the stationary phase approximation the GH shift will be positive. 
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Figure 5. (Color online) Lateral shift as a function of the angle of incidence for 
d/X = 0.7 and different widths of the gaussian incident wave. 


4. Transverse structure of the reflected beam 

To verify the prediction from the zero-pole model and the stationary phase 
approximation that the sign of the GH shift near the excitation of a backward surface 
wave depends critically on the thickness of the coupling layer, we compare the results 
of the previous section with those obtained from the rigorous calculation of the spatial 
structure of the reflected beam (equation ([5])). The integrals in equation (|5]) and ([6]) 
have been numerically calculated. 



Figure 6. (Color online) Lateral shift as a function of the angle of incidence for 
d/X = 0.9 and different widths of the gaussian incident wave. 

Figure |5] (for d/X = 0.7, a value lower than the critical value) and[6](for d/X = 0.9, 
a value higher than the critical value) show curves of the normalized hrst moment 
of the reflected electric held given by equation ([6]), as functions of the angle of incidence 
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9 and for different values of the beam width a. In all the cases, the sign of Ad) is in 
complete agreement with the zero-pole model and the stationary phase approximation. 
For both values of d/X, the larger the width of the incident beam, the greater the peak 
value of |Ad)|. Due to the shift of the SPP resonance with the coupling layer thickness, 
peak values occur for different angles of incidence {6 ^ 48° for d/X = 0.7 and 6 ~ 47.9° 
for d/X = 0.9). Note that large peak values for |Ad)| are obtained in both cases. When 
d/X = 0.7, the peak values of Ad)/A are —61.18, —89.19 and —97.02 for a/X = 50,100 
and 200 respectively, whereas when d/X = 0.9, the peak values of Ad)/A are 8.67, 41.18 
and 80.53 for a/X = 50, 100 and 200 respectively. 

Under conditions as those considered in these examples, such as total absorptions 
and resonant excitation of SPPs, the transverse structure of the reflected beam can 
have a complicated form. Particularly, the spatial distribution in the beam can be 
asymmetric, so that the value of the shift dehned by equation ([T]) may not coincide 
with the value obtained from the hrst moment of the spatial structure of the reflected 
beam. In these cases, and following a procedure similar to that presented in a paper 
by Tamir and Bertoni [ 8 ], it is possible to obtain a description which i) accounts for 
the deformation of the beam, ii) agrees quantitatively with the numerical calculations, 
and iii) takes advantage of the phenomenological model described by equation ([7]). The 
advantage of using the zero-pole model to represent the reflection coefficient is that the 
integrals in equation ([5]) can be evaluated explicitly in terms of the complementary error 
function of complex argument | 8 ]. In this description, the electric field of the reflected 
beam at the interface y = 0 results 


Er{x, 0) = C Ei{x, 0 ) 1-i TT'^' ni{zp - zq)— exp( 7 ^) erfc(- 7 ) 


( 8 ) 


where ni = y/eijxi and erfc is the complementary error function evaluated at the complex 
argument 7 

X .ax , s 

7 =- V i—'Kni{z - Zp). (9) 

Thus, the reflected beam is obtained as a product of two terms: one with a Gaussian 
form, which is multiplied by another function so that the reflected field is no longer 
Gaussian. To gain a deeper understanding of the distortions of the reflected beam, we 
have calculated the transverse distribution of the reflected electric field as a function of 
the a:-coordinate. As shown in Figures [7] and El the results obtained using equation ([5]) 
(solid lines) are in very good agreement with those obtained using equation ([ 8 ]) (circles). 

Figure [7] shows \Er{x,0)\ (equation ([5])) for d/X = 0.7, 6 = 48.01° (the angle of 
incidence for which the peaks in Ad) occur, see figure ED and for different values of 
the beam width a. The spatial distribution of the incident electric field (dashed line) 
is given as a reference. In these examples we observe that the spatial distributions 
are always shifted toward negative values of x, in total agreement with the values of 
A and Ad) obtained for a structure with the coupling layer thickness lower than the 
critical value. However, the profiles of the reflected fields are distorted with respect to 
the Gaussian profile of the incident beam, and the more severe distortion is obtained 
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Figure 7. (Color online) Profiles of \Er{x, 0)| calculated with equation ([5]) (solid line), 
\Er{x, 0)1 calculated with equation ([8]) (circles) and the Gaussian incident field (dashed 
line) for d/X = 0.7, 9 = 48.01° and values of a/X considered in figure O 
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for the smallest beam width. The prohles exhibit a splitting of the reflected beam. As 
explained in [3T], the major peak is produced by the spectral components of the incident 
beam which are coupled to the excited SPP, while the minor peak is produced by the 
spectral components of the incident beam which are not coupled with SPPs. 

For a structure with the coupling layer thickness higher than the critical value, the 
spatial distribution of reflected electric held is expected to be shifted toward positive 
values of x. The rigorous calculation of \Er{x,y = 0)| conhrms this expectation, as 
shown in hgure [HI analogous to hgure [71 but for d/\ = 0.9 and 9 = 47.86° (the angle 
of incidence for which the larger values for occur, see hgure [6]). Except for the 
sign of the lateral shift, the features exhibited by the prohles in hgure [HI are completely 
analogous to those in hgure [3 As in the previous example, a double peak structure in 
the rehected held prohle is present and the more severe distortions occur for the smallest 
beam width. 

5. Conclusions 

We have presented an exhaustive study of the GH lateral shift in Otto-ATR systems with 
negative index metamaterials near conditions of resonant coupling between the incident 
held and backward surface waves. Although the GH beam shift at a single flat interface 
can be either positive or negative depending on the type of surface waves (forward or 
backward) excited by the incoming beam, the erroneous belief that the excitation of 
a backward surface wave always results in a negative GH shift frequently appears in 
the literature, particularly when multilayers with negative index materials are involved. 
To clearly evidence that the excitation of backward surface waves in multilayers with 
negative index metamaterials can lead to both negative and positive (and not exclusively 
negative) GH shifts, we have adapted to the Otto conhguration the zero-pole model 
presented in [2l] for the Kretschmann conhguration. This model predicts the existence 
of a critical thickness of the coupling layer: above this critical value the lateral shift is 
positive whereas below this critical value the lateral shift is negative. We have presented 
rigorous calculations of the transverse distribution of the rehected electric held near the 
resonant excitation of the backward surface plasmon polariton. These calculations show 
the existence of large (negative or positive) Goos-Hanchen shifts accompanied by a 
splitting of the rehected beam and conhrm the predictions of the zero-pole model. 
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Figure 8. (Color online) Profiles of \Er{x, 0)| calculated with equation ©(solid line), 
\Er{x, 0)1 calculated with equation © (circles) and the Gaussian incident field (dashed 
line) for d/X = 0.9, 9 = 47.86° and values of a/X considered in figure El 
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